On Constructing the Asymptotic Solutions for 
Phase Transitions in a Slender Cylinder 
Composed of a Compressible Hyperelastic 
Material with Clamped End Conditions 

Hui-Hui Dai^, Jiong Wang^ and Zhen Chen^ 
^Department of Mathematics and Liu Bie Ju Centre for Mathematical Sciences, 
City University of Hong Kong, 83 Tat Chee Avenue, Kowloon Tong, Hong Kong 

Email: mahhdai@cityu.edu. hk 
^Department of Mathematics, City University of Hong Kong, 83 Tat Chee Avenue, 

Kowloon Tong, Hong Kong 
Email : j iongwang2@st udent . cityu . edu . hk 
^Department of Civil and Environmental Engineering, 
University of Missouri-Columbia,Columbia, MO 65211-2200, USA 

Abstract 

In this paper, we study phase transitions in a slender circular cylinder 
composed of a compressible hyperelastic material with a non-convex strain 
energy function. We aim to construct the asymptotic solutions based on an 
axisymmetrical three-dimensional setting and use the results to describe the 
key features (in particular, instability phenomena) observed in the experi- 
ments by others. The difficult problem of the solution bifurcations of the 
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governing nonlinear partial differential equations (PDE's) is solved through a 
novel approach. By using a methodology involving coupled series-asymptotic 
expansions, we derive the normal form equation of the original complicated 
system of nonlinear PDE's. By writing the normal form equation into a first- 
order dynamical system and with a phase-plane analysis, we manage to deduce 
the global bifurcation properties and to solve the boundary- value problem an- 
alytically. The asymptotic solutions (including post-bifurcation solutions) in 
terms of integrals are obtained. The engineering stress-strain curve plotted 
from the asymptotic solutions can capture the key features of the curve mea- 
sured in a few experiments (e.g., the stress drop, the stress plateau, and 
the small stress valley). It appears that the asymptotic solutions obtained 
shed certain light on the instability phenomena associated with phase tran- 
sitions in a cylinder, in particular the role played by the radius-length ratio. 
Also, an important feature of this work is that we consider the clamped end 
conditions, which are more practical but rarely used in literature for phase 
transition problems. 

Key words: Phase transformation; Hyperelastic material; Asymptotic analysis; 
Cylinder; Bifurcations of PDE's 

1 Introduction 

Applications of phase-transforming materials such as shape memory alloys (SMAs) 
and shape memory polymers are very broad. For example, they have been used to 
make satellite dampers, golf club heads and snake-like robots and so on. In particu- 
lar, these materials have been used to design many minimal surgery devices (Pelton 
et al. 1997). A deep and thorough understanding of the behaviour of this type of 
phase-transforming material is essential in the manufacturing and designing of these 
devices. 
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Systematic experiments on uniaxial extensions of superelastic NiTi alloy (a kind 
of shape memory alloys and also one kind of phase-transforming materials) wires and 
strips (Shaw & Kyriakides 1995, 1997, 1998) showed the measured engineering stress- 
strain curves have the key features: The nucleation stress occurs at a local maximum 
which is significantly larger than the Maxwell stress; following the nucleation stress 
there is a sharp stress drop; and afterwards there is a stress plateau. These features 
were also observed in experiments done by others (Sun et al. 2000, Tse & Sun 2000, 
Favier et al. 2001 and Li & Sun 2002). 

Theoretically, solid-solid phase transitions have also been studied for a long 
time in the context of both continuum and lattice theories. The seminal work 
of Ericksen (1975), which considered a continuum one-dimensional stress problem, 
made clear that for a non-convex strain energy function the solution with two phases 
can arise and there are multiple solutions. Based on the lattice model for a two- 
phase martensitic material, it is possible to deduce in the related continuum model 
the strain-energy function has double wells; see Ball & James (1992). In general, it 
is now understood that in the continuum scale for a material whose strain-energy 
function is non-convex phase transitions can take place (see also, Abeyaratne et 
al. 2001). With a proper choice of a strain-energy function, phase transitions can 
be modeled through a continuum theory. Still, to justify this point of view, it is 
desirable (probably necessary) to compare the analytical solutions based on this 
type of energy functions with experimental results. But, the difficulty is that we 
lack mathematical theories for mixed type equations that typically arise from non- 
convex energy functions. Analytical solutions for boundary- value problems are very 
few. In the classical paper of Ericksen (1975), analytical solutions were constructed 
for a static problem based on a pure one-dimensional stress model, which neglects 
the effect from other dimensions. However, a one-dimensional model appears to be 
not sophisticated enough to capture some key features observed in experiments. We 
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give the explanations below. 

In the experiments (Shaw & Kyriakides 1995, 1997; Li & Sun 2002), it was 
observed that after the nucleation stress was reached the nucleation process began 
and accompanying with it there was a radial contraction (necking). Also, after the 
two-phase state was formed, the deformation was inhomogeneous, one part being 
thin and one part being thick. As pointed out by Chang et al. (2006) that the axial 
extent of the transformation front is of the order of the radius. Thus, to model the 
nucleation process and the inhomogeneous deformation of different thicknesses, one 
should consider the radial deformation. 

Theoretically, if one treats the problem as a one-dimensional stress one, there 
is a discontinuous interface between two phases with two different strain values. 
The thicknesses of the cylinder at two phases are different, thus the shear strain Uz 
{U is the radial displacement and Z is the axial coordinate) is nonzero (infinite at 
the phase boundary). The traction-free boundary conditions require the two stress 
components S/j^ (depending on f/| up to the second order) and S/j^ (depending on 
Uz) to be zero at the lateral surface, which cannot be satisfied. Thus, it appears 
that to model phase transitions by a one-dimensional stress approximation violates 
the traction-free boundary conditions. 

In the papers of Dai & Cai (2006) and Cai & Dai (2006), phase transitions in a 
slender cylinder composed of a special incompressible elastic material was considered. 
A novel series-asymptotic approach is utilized to reduce the field equations. A proper 
asymptotic model equation is derived, which takes into account the influences of the 
radial deformation and traction-free boundary conditions. The solutions for two 
boundary-value problems are obtained, and they could capture the key features 
observed in experiments. 

In this paper, we study the phase transitions in a slender cylinder composed 
of a general compressible elastic material due to tension/extension. Also, different 
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from the previous studies on incompressible materials, here we further consider the 
couphng effect of the material nonlinearity and geometrical size. Another important 
new feature is that we consider the more practical clamped end conditions instead 
of the natural boundary conditions used in Cai & Dai (2006). 

We consider the problem in a three-dimensional setting, different from the one- 
dimensional stress problem studied by Ericksen (1975). However, the strain-energy 
function is assumed to have the same property as that in Ericksen (1975), i.e., for a 
one-dimensional stress problem the stress-strain curve has a peak- valley combination 
(cf. figure 1). We aim at constructing the asymptotic solutions and using them to 
explain the experimental results. 

Mathematically, to deduce the analytical solutions for the present problem is 
a very difficult task. One need to deal with coupled nonlinear partial differential 
equations (PDE's) together with complicated boundary conditions. Further, the ex- 
istence of multiple solutions (corresponding to the instability phenomena (e.g., stress 
drop) observed in experiments) makes the problem even harder to solve. Here, the 
analysis is carried out by a novel methodology developed earlier (Dai & Huo 2002; 
Dai & Fan 2004; Dai & Cai 2006), which is capable of treating the global bifur- 
cation problem of nonlinear PDE's and obtaining the post-bifurcation solutions. 
We construct the solutions and extract from them important information on the 
deformation configurations, the nucleation stress, the instability phenomena and 
the transformation front. Comparisons with experimental results are made, which 
show that the asymptotic solutions can capture the key features of the experimental 
engineering stress-strain curves and the instability phenomena as observed in ex- 
periments. The qualitative agreements give supporting evidence that a non-convex 
strain energy function can be used to describe solid-solid phase transitions and the 
instability phenomena in phase transformations are mainly due to the non-convexity 
of the strain energy function. 
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This paper is arranged as follows. In section 2, we formulate the field equations 
by treating the slender cylinder as a three-dimensional object. In section 3, we carry 
out a non-dimensionalization process to extract the important small variable and two 
small parameters which characterize this problem. Then we derive the normal form 
equation of the original governing nonlinear PDE's in section 4, through series and 
asymptotic expansions. In section 5, we show that the Euler-Lagrange equation can 
also lead to the same normal form equation, which justifies our method in deriving 
this equation. In section 6, we propose the clamped boundary conditions. In section 
7, we construct the asymptotic solutions for both a force-controlled problem and a 
displacement-controlled problem. We also use the solutions obtained to explain some 
experimental results. Finally, some conclusions are drawn. 

2 Three-dimensional field equations 

We consider the axisymmetric deformations of a slender elastic cylinder sub- 
ject to a static axial force at two ends. The lateral surface is traction-free and the 
end conditions will be considered later. In an undeformed state, the radius of the 
cylinder is a and the total length is 21. It is assumed that 6 = a/l << 1. We take 
the cylindrical polar coordinate system and denote {R, 6, Z) and (r, 6, z) the coordi- 
nates of a material point of the cylinder in the reference and current configurations, 
respectively. The finite radial and axial displacements can be written as 



We introduce the orthonormal bases associated with the cylindrical coordinates and 
denote these by E/j, Ee, E^ and e^, e^, in the reference and current configura- 





(2.1) 



6 



tions, respectively. Then the deformation gradient tensor F is given by 

F = {l + UR)er®ER + Uzer(8)Ez + il + ^)ee(8)Ee + WRe,0ER + {l + Wz)e,^Ez. 

(2.2) 

For an hyperelastic material, the strain energy function $ is a function of the three 
invariants Ji, I2 and of the left Cauchy-Green strain tensor B = FF"^; that is, 
$ = I21 h)- We suppose that $ is non-convex in a pure one-dimensional stress 
problem such that phase transition can take place. The nominal stress tensor S is 
given by 

^"9F' ''~dF\i ^^'^^ 

If the strains are small, it is possible to expand the nominal stress components 
in term of the strains up to any order. The formula containing terms up to the third 
order material nonlinearity is (cf. Fu & Ogden 1999) 



^ji = (^]ilkVkl + ^a^lknmVklVmn + ^a^lknmgpVklVmnVpq + 0{\r]st\^), (2.4) 



where rjij is the components of the tensor F — I, a^m., ajuknm o^jUknmqp ^'^^ 
incremental elastic moduli, which can be calculated once a specific form of the 
strain-energy function is given, and the formulas can be found in Appendix A. From 
the formula (2.4), we could obtain the nominal stress components Sj,. For example, 

U 1 f/2 „ UWz 9 UUr 

Sz. =i2^ + iiWz + i2UR + -(r/2-^ + r/4t/| + 2^/2-^ + + 2^/3-^ 



+ '2v2WzUr + ti^U'r + 2t^jUzWr + ii^W'^) + -(^2^ + 30, 



+ 3^3^^ + ^e^ulWz + 3^2^ + e,wl + 3^4^^ + ^e^ulUR 



R R R 

+ 69uUzWzWr + QOuUrUzWr + S^g— ^ + 3e,W^Wz + se^URW^), 

ri 

(2.5) 



where ^j, rjj and 6k are some elastic moduli, whose formulas are given in Appendix A. 
The other non-zero stress components can also be obtained but we omit their lengthy 
expressions for brevity. Owing to the complexity of calculations, we shall only work 
up to the third-order material nonlinearity. In the experiments, the maximum strain 
is less than 10%, and such an approximation is sufficiently accurate. 

S satisfies the following field equations: 

Div(S) = 0, 

which yields the following equations 

d^Zz dY^Rz Yrz _ 
dZ dR ~W ~ 

dT^Rr dT^Zr ^Rr — ^09 

dR dZ R 

We consider the case that the lateral surface of the cylinder is traction-free 
we have the boundary conditions 

S_Rr|_R=a = 0, Sfc|iJ=a = 0. (2.9) 

We shall derive the asymptotic solutions of (2.7) and (2.8) under (2.9) and some 
end conditions for a given external force. 

Equations (2.7) and (2.8) together with (2.4) provide the governing equations 
for two unknowns U and W. The formal two are very complicated nonlinear partial 
differential equations (PDE's) and the boundary conditions (2.9) are also compli- 
cated nonlinear relations (cf. (3.3)-(3.6)). To describe the instability phenomena 
during phase transitions, one need to study the bifurcations of this complicated 
system of nonlinear PDE's. Here, we shall adopt an alternative approach involving 
coupled series-asymptotic expansions to tackle this bifurcation problem. A similar 
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2.6) 



(2.7) 

(2.8) 
Thus, 



methodology has been developed to study nonlinear waves and phase transitions in 
incompressible materials (see Dai & Huo 2002, Dai & Fan 2004, Dai & Cai 2006, Cai 
& Dai 2006). First, we shall nondimensionalize this system to identify the relevant 
small variable and small parameters. 

3 Non-dimensional equations 

We first introduce a very important transformation 

U = vR, s = i^^ (3.1) 

then do the following scalings: 

s = l^s, Z = lz, W = hw, v = jv, e = y, (3.2) 

where / is the length of the cylinder, /i is a characteristic axial displacement, and 
e is regarded to be a small parameter (equivalent to a small engineering strain). 
Substituting (3.1) and (3.2) into (2.7) and (2.8), we obtain 

2(6 + ^3)v, + iiw,, + 4^3^. + ^[2(6 + ^:i)vsz + 4^3^;.^] + ■ ■ ■ = 0, (3.3) 

i^v,, + SiiVs + 2(^2 + iz)wsz + sA^iVss + • • ■ = 0. (3.4) 

Here and hereafter,we have dropped the tilde for convenience. The full forms of (3.3) 
and (3.4) are very lengthy and can be found in Appendix B, here we just present 
the first few terms. Substituting (3.1) and (3.2) into the traction-free boundary 
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conditions (2.9), we obtain 

13 1 

(6 + ^2)v + ^2Wz + s2^iVs + 4(2'^! + 2^12^ + (^72 + ni)vw-^ + -ri2wl 

+ si^V^vl + (2?7i + 2t]2)vVs + 2r]2WzVs + 2r]YV^Ws + 2r]4wl) + s^2t]iv'^] 

+ + ^^2 + ^^3)^' + (^^2 + ^9,yw. + il-Os + l9,)vwl + ]:92wl (3.5) 

boz Z Z Z Z 

+ S((^^5 + ld8)vvl + ^^7^^>. + (^1 + 202 + e^yVs + {202 + 2^4)t^«^.^^s 

+ Oswlvs + (2^12 + 2^i5)i;i;^u;s + 2euv^WsW^ + (2^5 + 2^8)t'u^s + 2e7W^wl) 



4 

+ s^(e5'y>. + (2^1 + 2^2)^^^^,^ + 292W,vl + A9i2VsV,Ws + AO^Vswl) + s^-eiv%=^ = 0, 

^3V^ + 2^3Ws + e[(?77 + r]s)vv;, + rjjv^^W:, + (2774 + 2ri^)vWs + 2?74W^w, 
+ s{2r]jv^Vs + 4?74t;,Ws)] + e^[(i6'i2 + ^6*13 + 6ir,)vh, + (6^14 + 6i5)vv^w, 

+ \o^2V,wl + (^5 + ^6 + 2^8)^^'w^s + (2^7 + 2^8)^w.«^. + O^wlws (3.6) 
+ s(^^i6t'^ + (2^12 + 2^15)1^^1;, + 2euv,w,Vs + Oi-jvlws + (4^5 + Ae^)vVsWs 



4 

+ 4^7w;,t;,u', + 2eiQV,wl + -^gUi.^) + s^{2ei2Vzv1 + 4^5W>,)]|,=^ = 0, 

where v = a^/t^ (square of the diameter-length ratio) is a small parameter for a 
slender cylinder. We note that in the above equations the dependence on R and a /I 
is entirely through s = R'^/P and u = a? /f', respectively. 

Equations (3.3)-(3.6) comprise a new system of complicated nonlinear PDE's 
with complicated boundary conditions, which is still very difficult to solve directly. 
However, it is characterized by a small variable s (0 < s < z/) and two small 
parameters e and z/, which permit us to use expansion methods to proceed further. 
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4 Coupled series-asymptotic expansions 

From equations (3.3)-(3.6), we can see that the two unknowns w and v depend 
on the variable z and the small variable s and the small parameters e and that is 

w = w{z, s;e,v), V = v{z,s]e,h'). (4.1) 

To go further, we seek series expansions in terms of the small variable s: 

v{z, s; e, v) = Vo{z; e, u) + sVi{z; e, u) + s'^V2{z; e,u) -\ , (4.2) 

w{z, s; e, u) = Wo{z; e, u) + sWi{z; e, u) + s^W^^z; e, i/) + ■ ■ ■ . (4.3) 

Substituting (4.2) and (4.3) into (3.3), the left hand becomes a series in s and all 
the coefficients of s" (n=0, 1, 2, ■ ■ ■ ) should vanish. Equating the coefficients of s° 
and to be zero yields that 

4^41^1 + (26 + 26)^. + ixW^zz + e((4r/4 + 4r/5)K)W^i 

+ (2r/2 + 2r/3 + l^i + 27/8)K)K). + ^'^^W^W^, + (2r/2 + I^tW^zW^z (4-4) 
+ 1^2y^W^zz + r/iW^o.W^o..) + e'i^i(^, W^o, W^i) = 0, 

16e3W^2 + (46 + 46)Vi. + 6W^i.. + e((24r/4 + 81/5)^1 W^i 

+ (16r/4 + \^rii)V^W2 + (4r/2 + 47^3 + 12r/7 + 4r/8)Vi^. 

+ (4r/2 + 4r/3 + 4r/7 + 4r/8)^^i. + 16r/4W^2W^o. + (47?2 + 4r/7)V^i,W^oz 

(4.5) 

+ ll^^WxWxz + (2r/2 + 6r/7)K).l^i. + 27/7^^1 + ^4^).^).. 
+ e2i72(K,,^i,iyo,H^i,W^2) = 0. 
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Similarly, substituting equations (4.2) and (4.3) into (3.4) and equating the coeffi- 
cient of s° to be zero yields that 

S^iVi + (2^2 + 2^3)1^1. + ^sVo.. + e((8r/i + 87/2)^0^1 
+ (6r/4 - 2r],)W^ + SrjrWiVo, + (^r/4 + ^V5)Vol 

+ 8r]2ViWo, + (2r/2 + 2r/3 + 2r/7 + 2r]s)VoWi, + (2r/2 + 27/7)1^0.^^1. (4-6) 

+ e'H,{Vo,V^,Wo,W,)=0. 

The expressions of Hi, if 2 and if 3 are very lengthy, which are omitted for brevity. 
Substituting (4.2) and (4.3) into the traction-free boundary conditions (3.5) and 
(3.6), we obtain 

(^1 + ^2)Vo + 6W^o. + i^mi + 6)^1 + 6^^!,) + e((^77i + |r/2)K,2 

^ ^ (4.7) 

+ {V2 + VsWoWo, + -r/2W^o'J + e'H,{Vo, Wo) + euHM, V^, Wo, W^) = 0, 

2^3^^! + ^sVoz + i^(46W^2 + ^3Vi,) + e((2r/4 + 2r]5)VoWi + (r/7 + r/s)^,^. 

+ 2r]^WiWo, + 777K).l^o.) + e'H,{Vo, Wo, W^i) + ei.ii7(\/o, V,, Wo, W,, W^) = 0, 

(4.8) 

where the lengthy expressions for H^-H-j are omitted, and interested readers can 
find the expressions for iij (i = 1, 7) in Appendix C. In obtaining the above 
equations, we have neglected terms higher than 0{ev,e^). 

Now the governing equations (3.3)-(3.4) are changed into a one-dimensional 
system of differential equations (4.4)-(4.8) for the unknowns Wo, Wi, W2, Vo and 
Vi. By using a regular perturbation method, we can express Wi by Wq and Vq from 
(4.4), then substitute Wi into (4.6), we can also express Vi by Wo and Vq. Similarly, 
we substitute Wi and Vi into (4.5) to obtain the expression of 14^2 in terms of Wo 
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and Vq through the regular perturbation method. The results are given below: 

^1 + Sr)^ozz + (4 + -^Wozzz + eianVol + aUVoVo,, 

44i 84i43 16 1643 

+ aiaWozVozz + OmVo^W^ozz + ^i^Wq,^ + ax^V^WQ^zz + o-x-jWqz^qzzz) 

+ anWl^V^zz + 023^)^0^1^02^ + a^AV^zW^zW^zz + a^^y^Wl^^ 
+ aseWo^W^o^^ + a27'K)^W^ozzz + oasKjW^ozW^ozzz + a^^Wl^W^zz-^, 

^ ^ 16^1 32^1 32^1^1 32^3 32^1^3^ ^ 64 64^| 

^2 

^ «772 ~ QOf y^Ozzzz + e('^3oK)z^02z + 0,3iVozzWqzz + d'^'l^O^Ozzz 
0443 d^43 

+ a33^0zV^)^^^ + 034^0^1^0222 + 035^^022^0222 + ^36^^02222 + 037^0^1^02222) 
+ e^(a38Vo^^ + OggVoK)^ 14)22 + a4oV'o2W^02K)^2 + 041^^2^022 + 042^)^)221^022 
+ 043^^02^)22 1^022 + a^^Voz^lzz + (^^^^Izz + 046K?K)222 + 047^^^021^222 

+ 048PV0^2K)222 + 0491^01^02^^0222 + 050 ^21^02^0222 + O5I VJ) 1^022 1^0222 

+ 0521^02 1^022 1^0222 + O53 1^)^1^02222 + 0541^1^02^^02222 + 055^^02^02222), 

(4.11) 

where Oj [i = 1, 2, ■ ■ ■ , 55) are constants related to material constants, whose ex- 
pressions are omitted for brevity. Substituting Wi, V\ and W2 into (4.7) and (4.8) 
and omitting the higher order terms yield the following two equations with only two 
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unknowns Wq and Vq: 

(ei + + 6W^o. + e{bX + hVoWo, + hW^,) + e^hV^^ + hV^Wo. 

+ hV^Wl + hWl) + //(fesK).. + fegW^o...) + ei/(6ioV;,', + 611^)^0.. (4-12) 

+ &i2Vro^K)^^ + bM^Wo,, + ftuVTo^^ + ^isVoVTo^^^ + hiQWoMozzz) = 0, 

and 

+ 6^(^21 K^V^O. + b22VoVozWoz + b2zVo,W^^ + 624Kf H^O.. + b25VQWo,Wo,, 

+ &26W^o,W^o,^) + T^{b27VQ,,, + 628^^0,,,,) + ez/(629M).V^o.. + b^ioVozzWo,^ (4-13) 
+ ^siVbVoz.. + b'i2Wo^Vo^^^ + 633^02^^02.. + b^iWozzWo^^^ + 635^0^^0222. 
+ 636^^0.^^02222) = 0, 

where 6j (i = 1, 2, ■ • • , 36) are also constants related to material constants, whose 
expressions are omitted for brevity. (Note: Interested readers can contact the corre- 
sponding author for the formulas for Oj and 6j.) By using the above two equations, 
we can express Vq in terms of Wq and its derivatives with respect to z. From (4.12), 
we have 

^0 =^^[6W^02 + e{b^V^ + b2VoW^, + 63^^02) + e'(&4^o' + hV^W^z 
u + « 

+ b^V^Wl + hWl) + //(fegVo.. + fegW^o...) + eu{bi^Vl + ^^"^^^ 
+ &i2Wo2K)2^ + 6i3K).W^o.. + &i4Wf;2^ + fei5K)W^o... + ftieW^o.W^o...)]- 

From (4.12) and (4.13), we also have 

^0 = ^^[6W^o. + e(6iK)' + ^2Wo. + Wo.) + K&8K).. + Wo...)] + --- , (4.15) 
4i + 42 
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and 



Substituting Vq = -^^^oz + ■■■ into (4.15) and (4.16), then substituting (4.16) 
into uip^Vozz + bgWozzz) of (4.14) and (4.15) into the other part of (4.14), we obtain 



6+« 16(4l + ?2) 

+ ei/(a3iyo2^^ + aiWozWozzz), 

where 



(4.17) 



2(^1 + 6) 



«2 , 05 [^1^^^2 + ^V2V3 - ^102) + 36^2(2^/2' + 47/2773 - 4r/| 



3 

1 

6(6 + 6) 

+ Vi{V2 + 3r/3) - 36^2 + 36^3) + ^2(^1 + 6^3 - 6^4) + ^li-Svl 
- 12r/ir/2 - 12r/2 + 9r/ir/3 + 18r/2r/3 - 6r/| + 6^1 - 66^2 + 126^3 



- 96^4) - eM^Vm + 6vl + 3r/2^3 + Gvl + 46^2 - 36^3 - 36^4)], 

«3 =77777— T^l^^^t + 4^2(4^1 - 3r/2) - 6'6(236 + r/i + 13r/2 - 267/3) 
646(6 +6) 

+ 46^2(3^/1 - 10772 + 4773) + e?(106 - 3771 - 7772 + 6773)], 
«4 =77^^77 ^.^ [3^1 + 26'(-m + ^2) + 6'e2(-516 + IO772 - 6773) 

166(6 -6) (6 + 6) 
+ e^6(-216 - 5771 + 772 - 4773) + 6'^2(-306 + 7771 + 21772 - 4773) 
+ 6'(36 -Vi + 5r72 + 2773) - 6el(7r7i + 157/2 + 4773)]. 
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By substituting (4.16) and (4.17) into (4.13), we obtain 

4 (4.18) 

+ e{2D^WMzz + y{^D^WQ,,Wozzz + ID^WozWozzzz)) = 0, 

where 



+ 36^2(^1-27/2 + 27/3)], 



(4.19) 



- 86^2) + ^2(3^1 - 8^2 + 18^3 - 12^4) + 2^f^j{-6riiri2 - Urij + 12ri2m 

- 127/3 ^ 5^^^^ _ ^Q^^Q^ ^ g^^^^ ^ g^^^^) ^ 2^1^2(127/17/3 + 247/27/3 

- 247/^ + 56^1 - 286^2 + 186^3 + 66^4) - ^2(67/^ + 24r]l + 247/ir/2 

- 247/17/3 - 487/27^3 + 24r/| - T^^i + 40^1^2 - 426^3 + 12^1^4)], 



(4.20) 

^3 =^77 AtT—TYsI^^i + ^1(4^2 - r/i - 3r/2) + 6e2(-46 - 7r/i + 3r/2) 

+ 6^2(^2 - ^3) - 2e?6(e2 + 27/1 - 37/2 - 37/3)]. 

(4.21) 

Integrating (4.18) once, we obtain 



eWoz-\euWo,,, + e^D,Wl + e^D2W^, + e'u{D3Wl^ + 2D3Wo,Wo,,,) = A (4.22) 



where A is the integration constant. It is important to find the physical meaning 
of A, since to capture the instability phenomena observed in the experiments, one 
needs to study the global bifurcation as the physical parameters vary. For that 
purpose, we consider the resultant force T acting on the material cross-section that 
is planar and perpendicular to the cylinder axis in the reference configuration, and 
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the formula is 

/'27r ra 

T= / Ez.RdRdO. (4.23) 
Jo Jo 

By using (4.2), (4.3) and the expressions of Vq, V\ and W\ in terms of Wq in (2.5), 
it is possible to express in terms of Wq. Then, carrying out the integration in 
(4.23), we find that 

(4.24) 

where E = Young's modulus. Comparing (4.22) and (4.24), we have 

A = If we retain the original dimensional variable and let V = Wqz = ^-Wqz-, 

we have 

V + D^V^ + D^V^ + a\~Vzz + D^Vi + 2D3^^zz) = 7, (4-25) 

where 7 = ^^^^ is the dimensionless engineering stress. Since (4.25) is derived 
from the three-dimensional field equations, once its solution is found, the three- 
dimensional strain and stress fields can also be found. Also, it contains all the 
required terms to yield the leading-order behavior of the original system. Therefore, 
we refer (4.25) as the normal form equation of the system of nonlinear PDE's (3.3) 
and (3.4) together with the boundary conditions (3.5) and (3.6) under a given axial 
resultant. 

Remark: Although the final form of the normal form equation is one dimensional, 
it is different from other one-dimensional phenomenological models in literature (see 
Coleman 1983, Tong et al. 2001, Shaw 2002, Chang et al. 2006) in the sense that it 
is derived from the three-dimensional field equations in a mathematically consistent 
manner and all three-dimensional quantities can be calculated once its solutions 
are found. Thus, we can provide the three-dimensional current configuration in a 
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two-phase state, which cannot be deduced in the existing one-dimensional models 
in literature. This equation contains a higher-order derivative term which plays 
the role of regularization. Trunskinovsky (1982, 1985) pioneered the idea of the 
regularization augmentation for solid-solid phase transitions which involves adding 
terms (like a strain gradient) into the usual constitutive stress-strain relation (for 
the gradient approach, see also Aifantis & Serrin 1983 and Triantafyllidis & Aifantis 
1986). Different from these gradient theories, the gradient term (which represents 
the influence of the radial deformation and traction-free boundary conditions) in 
our equation is derived and its coefficient is explicitly given. 



5 The Euler- Lagrange equation 

The normal form equation (4.25) is derived in the previous section based on the 
equilibrium equations (2.7) and (2.8). Now, we consider the energy and show that 
the same equation can be derived through the variational principle. The calculations 
are tremendously complicated in the present method (we use Mathematica for all 
the symbolic computations), and if the two derivations of the normal form equation 
agree with each other, it also supports the correctness of our calculations. 

We expand the strain energy function up to the fourth-order nonlinearity (which 
implies that the stress components are up to the third-order nonlinearity), and as a 
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result we obtain 

1 U IP 

+ 2^:,UzWn + esVT^) + -(ryi— + 3r/5-t/| + 3r/2-^W^z + ^VaUIWz 

+ 3r/2^t/^ + SV2WzUj, + r/i?7| + Gvs^UzWr + Qi^UzWzWr 
ri JrC 

+ GvjUzUrWr + 3r/5^1^R + ^ViWzW^ + Sti.UrW^) 

ri 

+ QOs^wl + ee.ulwl + 4^2^p^I + ^iW^I + QO^wlul 

+ m.^wlUR + ao^wIUr + 6e,^u'^ + Qe^ulul + 12^4^1^^^/^ 

+ 4^2^^| + ^e^WzUl + ^if/^ + 12^i3^f/zW^fi + 2A9^,UzWzUrWr 

+ 24^i5^?7ziyziyij + m^^UzWlWl + 24^i5^t/zf/i?W^i? + 4^i6[/|iy^ 
ri ri 

+ 12^i2t/zf/RW^z + QoJ^Wl + 6ei7f/|lV|, + 12^8^W^zW^^ + QO^WlWl 

K K 

+ uOs^UrW^ + m^WzURWl + Qe^ulwl + ao^^UzWI + ^gW^^)- 

(5.1) 

The strain energy per unit length is given by 

^RdRdQ. (5.2) 



By using the same manipulations in section 3 and section 4, we can get the average 
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strain energy over a cross section 



•qj = 



6 24 4^3 2^1 4^3 4^i43 



(5.3) 



4 2^1 4^3 4^1^' 4^3 166 



3 



Ozz 



+ (| + I - H + §:)yoWo.zz + (| - ^ + ^)W^o.W^o...) 

+ CeVoVozWozz + CrVo.Wo.Wozz + csVqW^,, + CgWo.Wl, + c^Vi^Wozzz 

+ CuVoWozWozzz + CuW^Wo,,,)], 

where q {i = 1, ■ ■ ■ , 12) are some constants, whose expressions can be found in 
Appendix A. 

By further using (4.12) and (4.13), we can reduce the above equation as 



+ -^D^e^W^, + eu{HiW^Ml, + H^WlW^zzz)] (5-4) 

--E[]-V^ + \d^^^ + -D^V^ - \a^VVzz + a\H,VVl + H^V^Vzz)]. 
16 4 o 



where 



^ e?(r/i + 2r/2) + 6i(3r/i + 2r/2 - 2r^^) + 2^^UVi - 2r/2 - 2r/3) 
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+ e2'(3r/i + 14(r/2 - r/3)) + 2e?6(36 - r/i + r/3) (^-S) 
+ eie2(-86 - 7r/i + llr/2 + 2r/3) - '2Ul{^i2 + 3r/i - 6r/2 + Tr/a)]. 
The total potential energy is then given by 

n= f mdZ-T I VdZ = TTa\[ ^>dZ-E f -fVdZ) 
J-i J -I J~i J~i 

=7ra'E f {-^V + + \d,V' + -D,V' - \a'VVzz (5.7) 
J-l 16 4 8 

+ a2(/7irV| + //2V^Vzz))f/^. 
Further by the variational principle, from the Euler-lagrange equation we obtain 



V + D^V'' + D^V^ + a\-U/zz + (2i/2 - + 2(2iJ2 - i^i) VV^zz] = 7- (5-8) 



It can be seen from (5.5), (5.6) and (4.21) that 2iJ2 — Hi = D3. So the equation 
(5.8) is exactly the equation (4.25). 

The expressions of and Q are themselves important. When there are multiple 
solutions, the smallest energy criterion can be used to judge the preferred solution 
(configuration). 



6 The clamped boundary conditions 

Now, we conduct a detailed analysis on the normal form equation with clamped 
end boundary conditions. We suppose that the strain energy function is non-convex 
in a homogeneous constant strain state such that phase transition can take place. 
This requires 

Di < 0, D2> 0, 3D2 <Dl< AD2. (6.1) 
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We rewrite equation (4.25) clS cL first-order system 



Vz = y 

_ V + Di^^ + D^V^ + D^aW - 7 ^^-^^ 
a2(l-2Z}3V) 

The critical points of tliis system are determined by y = and 

V + DiV^ + D2V^ = (6.3) 

In tliis paper, we always choose Di = —18, D2 = 100. We can see that with 
these material constants, the critical stress values 71, 72, 7m (cf. (6.4)) and the 
corresponding strain values (cf. Figure 1) are close to the experimental results (Tse 
& Sun 2000 and Li & Sun 2002). The 7 — curve corresponding to equation (6.3) 
is shown in Figure 1. 



Figure 1. The j — V curve corresponding to (6.3). 
The peak stress value 72, the valley stress value 71 and the Maxwell stress value 7„ 



are 



2D^ - 2(D^ - 3^2)^ - 9D1D2 
7i = — i ^ ^ — = 0.015929070313677, 

1^ ( J^n 



'2 



72 



^ — = 0.017670929686323, (6.4) 



27D 



2Df - 9D1D2 
f 

The natural boundary conditions Vz = at the two ends have been used by 



7^ = ^ = 0.016800000000000. 

27Do 
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many authors (e.g., Ericksen 1975; Tong et al. 2001; Cai & Dai 2006). However, 
such conditions are difficult to be reahzed in practice. In experiments, usually the 
two ends are clamped into the rigid bodies. To better describe the realistic situation, 
in this paper we will consider the case that the two ends of the cylinder are clamped 
into rigid constraints. In this case, we can propose the conditions according to the 
fact that there is no radial displacement for the point at the lateral surface of the 
ends, i.e., 

U\R=a = ^,at Z = -l,l. (6.5) 
By using (3.1), (3.2) and (4.2), neglecting the 0(i/^) terms, we obtain 

Vq + uVi = 0,at z = -/, /. (6.6) 

In deriving equation (4.18), Vq and Vi have been expressed in terms of Wq. Using 
such relations in equation (6.6) and with the help of equation (4.25), we obtain 

FiV + F2V^ + FsV^ = -f, at Z = -IJ. (6.7) 

where Fi, F2 and F3 are material constants. Notice that in deriving (6.7) we have 
neglected the terms higher than 0(e^, u). It's clear that if we choose different values 
of Fi, F2 and -F3, the properties of the boundary condition (6.7) may also have a 
lot of differences. Here we will only consider a simple case. We will choose Fi, F2 
and F3 such that the cubic equation (6.7) has only one unique positive root Vg when 
7 > 0, and this root is much smaller than the positive roots of equation (6.3). In 
this paper, we will choose Fi = 3, F2 = —20 and F3 = 100 to get the following 
graphic result: 
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Figure 2. Relationship between the ^ — V curves corresponding to (6.3) and (6.7). 
Thus the boundary conditions at the two ends are 



We note that Ve depends on the value of the engineering stress 7. 

Remark: Although the normal form equation is one-dimensional, it is derived 
from a three-dimensional setting, and as a result we can use the end conditions 
by considering the quantities in the radial dimension. If one directly introduces a 
one-dimensional model, such an option is not available. 

7 Asymptotic solutions 

We now consider the solutions for a given engineering stress 7 under the clamped 
boundary conditions (6.8). In the following, without loss of generality, we take 
/ = 1. Then in equation (4.25), a is equivalent to the diameter-length ratio. Due to 
symmetry, we only consider the part < Z < 1 and boundary conditions (6.8) can 
be replaced by 



V = Ve at Z 



I J. 



(6.8) 



Vz\z=o = 0, V\z=i = V,. 



(7.1) 
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7.1 Force-controlled problem 

We regard the engineering stress 7 as the bifurcation parameter. As 7 varies, 
there are seven types of phase planes, which are shown in Figure 3. 




Every trajectory in the phase planes of Figure 3 is a solution of the system 
(6.2). However, only those trajectories that satisfy equation (7.1) represent the 
physical solutions for the present problem. It can be seen that for a trajectory to 
be the physical solution a necessary and sufficient condition is that it contacts the 
\^-axis once (so that Vz\z=o = 0) and contacts the vertical line V = Ve once (so that 
V\z=i = Ve) and the Z-interval for this segment of the trajectory is exactly equal to 
1. 

To deduce the solution, we integrate (4.25) once to obtain 



H-jV + ^V' + + ^D^V + a\-^V^ + DsVV^) = 0, (7.2) 



where H is the integration constant. If one trajectory contact the V-axis at point 
{vq, 0), Vq must be a real root of the equation 



H--fV+lv^ + \DiV^ + ^02^ = (7.3) 



By using (7.2) and consider the case < Z < 1, we can get 



dV__ H - 7V + ^V^ + ID.V^ + ^D2V^ 



(7.4) 



with P = Here, for simplicity, we choose D3 = —5 such that P < for graphic 
results. Suppose that V\z=o = vq, we can obtain 



^ = -^J\ L._±Lr + M2.r^ + r4 ^^> V<V<v,^ (7.5) 



where 



= -(-7^0 + + ^D.v^o + |Awo)- 
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To satisfy the end boundary condition (7.1), we need to choose Vq such that 



Based on the phase planes in Figure 3 and equation (7.5) and (7.6), we can construct 
all the possible solutions for a given 7. 

Case (a) < 7 < 7^ (cf. Figure 3(a)-3(d)) 

In this case, from Figure 3(a)-3(d), we can see that only the trajectories located 
between the left saddle point Vi and the vertical line V = Ve (notice that Ve <^ Vi, 
cf. Figure 2) can be the possible solutions. 

For any given 7, and Vi can be calculated immediately from the two equations 
(6.3) and (6.7). By using the equation (7.6), we can determine the value of Vq. Once 
Vq is determined, through numerical integration, the corresponding solution can be 
obtained from (7.5). 

For any 7 satisfies < 7 < 7^, we found that there is only one solution (denoted 
by 5*1), which is very close to the saddle point Vi. 

Case (b) 7™ < 7 < 72 (cf. Figure 3(e)) 

In this case, for some given 7, there are multiple solutions. 

From figure 3(e), we can see that the trajectories located between the left saddle 
point Vi and the vertical line V = Ve can be the possible solutions. Through some 
calculations, we found that for any 7m < 7 < 72? there exist one solution (also 
denoted by Si), which is represented by a trajectory to the left of (also very close 
to) the saddle point Vi. 

We can see that the trajectories located between the homoclinic orbit and the 
right saddle point V3 can also be the possible solutions. Through some calcula- 
tions, we found that there is a critical engineering stress 7^ = 7p(a) which satisfies 
7m < 7p(o) < 72- When 7p(a) < 7 < 72, there exist another two solutions. One so- 




(7.6) 
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lution (denoted by 5*2) is represented by a trajectory outside (but very close) to the 
homoclinic orbit. The other solution (denoted by S3) is represented by a trajectory 
outside the second one and to the left of (very close to) the right saddle point V3. 

Thus the critical engineering stress 7p(a) can be considered as a bifurcation 
point. When 7m < 7 < 7p(o), there is only one solution. When 7p(a) < 7 < 
72, there are three solutions. We found that 7p(a) is a monotonically increasing 
function of the diameter-length ratio a, e.g., 7^(0.03) = 0.0168005239556, 7^(0.06) = 
0.0168907218423. 

Case (c) 7 = 72 (cf. Figure 3(f)) 

In this case, there are three solutions. The first solution (also denoted by Si) 
is represented by a trajectory to the left of (also very close to) the cusp point V2. 
The second solution (also denoted by ^2) is represented by a trajectory to the right 
of (also very close to) the cusp point V2. The third solution (also denoted by 5*3) is 
represented by a trajectory to the left of (also very close to) the saddle point V3. 

Case (d) 7 > 72 (cf. Figure 3(g)) 

In this case, there exist another critical engineering stress 7^ = 7q(a) which is a 
little bit larger than 72. 

When 72 < 7 < 7g(a), there are three solutions. The first two solutions (also 
denoted by Si and 5*2) are represented by two trajectories (very close to each other) 
located in the middle part of the interval {Ve, V3). The third solution (also denoted 
by 53) is represented by a trajectory to the left of (very close to) the saddle point 
V3. When 7 > 7q(a), there is only one solution (also denoted by S3), which is 
represented by a trajectory to the left of (very close to) the saddle point V3. 

Thus the critical engineering stress 7g(a) can be considered as another bi- 
furcation point. We found that 7g(a) is also a monotonically increasing func- 
tion of the diameter-length ratio a, e.g., 7^(0.03) = 0.017677609814, 7^(0.06) = 
0.017735345200. 
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For the radius a = 0.03, the axial strain distribution V{Z) corresponding to the 
above four cases are plotted in Figure 4. 




z 



(a) (b) 






id) 



Figure 4. The axial strain distribution curves for four different 7 values (a = 0.03): (a) 
< 7 = 0.010080 < 7p, (b) 7p < 7 = 0.017366 < 72, (0)7 = 72, 
(d) 7 = 0.020114 > 7g. 



From the above results, it can be seen that there are multiple solutions for 
7p(a) < 7 < 7g(a). We suppose that the preferred solution should be the one which 
has the smallest total potential energy value. To determine which solution is the 
most preferred one, we consider the total potential energy values of all the possible 
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solutions. From (5.7) and due to symmetric, we can get 



Jo 2 6 A 6 

+ a^H^VVl + H2V^Vzz))dZ, 



(7.7) 



Without loss of generality, we choose i^^i = 15 and H2 = 5 for graphic results. For 
the radius a = 0.03, we plot the total potential energies for all the solutions in Figure 
5. 

Q (27ra^E) 
-0.00018 

-0.0002 

-0.00022 

-0.00024 




Figure 5. Relationships between 7 and the potential energies (a = 0.03). 

From figure 5, it can be seen that there is a critical stress value 70 = 70(a) (e.g., 
7o(0.03) = 0.01694471757). For 7p < 7 < 70, the first solution 5*1 is the preferred 
solution, and for 7o < 7 < 7^, the third solution 6*3 is the preferred solution. Thus, 
in a loading process, after 7 > 70, the configuration of the cylinder may jump 
from solution Si towards solution 5*3, which is corresponding to the phase transition 
process. 

Remark: The existence of 70(0.) (> 7p(a) > 7m) is significant and it implies that 
under clamped end conditions the phase transition can only happen at a stress value 
larger than the Maxwell stress. 
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7.2 Displacement-controlled problem 

We now consider the case that the total elongation 

J VdZ = 2 VdZ = 2A (7.8) 

is given. Since we have taken / to be 1, A is actually the engineering strain. The 
governing equation is still equation (4.25), but now 7 is an unknown parameter. 

In the previous section, we have obtained all the solutions for a given 7. If for 
a given A we can find the corresponding value 7, then we can obtain the solutions 
for a displacement-controlled problem. 

We first plot the 7 — A curves corresponding to the solutions we obtained in 
the previous section for four different values of a in Figure 6. 




(a) 
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Figure 6. The engineering 7 — A curves for different values of a: (a) a = 0.02; (b) 
a = 0.03; (c) a = 0.04; (d) a = 0.06. The dashed hues are the curves shown in Figure 2. 

We can see that the 7 — A curves shown in Figure 6 have some important 
features. 

First, these curves obtained from the asymptotic solutions capture the main 
characteristics of the engineering stress-strain curves measured in a number of ex- 
periments (Shaw & Kyriakides 1995, Sun et al. 2000, Tse & Sun 2000, Favier et al. 
2001 and Li & Sun 2002), i.e., there is a stress peak, which is followed by a sharp 
stress drop, and then there is a stress plateau (cf. Figure 6(c)). 

Second, it can be seen that there is a snap-back when the diameter-length ratio 
a is small (cf. Figure 6(a) and 6(b)). As the diameter-length ratio become larger, 
this snap-back phenomenon disappear (cf. Figure 6(c) and 6(d)). 

Third, as the diameter-length ratio a increases, the whole 7 — A curve moves 
towards left, especially in the high-strain region. It is clear that this feature is due 
to the boundary effect. 

With these curves, we are ready to convert the solutions of a force-controlled 
problem into those of a displacement-controlled problem. As we pointed out that 
if the diameter-length ratio a is large enough, there is no snap-back on the corre- 
sponding 7 — A curve (cf. Figure 6(c) and 6(d)). Thus, in this case, for a given 
displacement A, we can only get one corresponding stress value 7 from the 7 — A 
curve. For this 7, there is a unique phase plane. There are maybe more than one 
trajectories in this phase plane which satisfy the end boundary conditions, however 
only one of them can give the required displacement value, which then represents 
the unique solution for the given A. 

As we pointed out before that the normal form equation (4.25) is derived from 
the three-dimensional field equations. With the solutions of the normal form equa- 
tion, we can recover the three-dimensional strain fields. Here, we take a = 0.04 as 
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an example. Corresponding to the four points A — D shown in Figure 7, 




Figure 7. The engineering 7 — A curve for a = 0.04. 

we draw the shapes of the cyhnder in Figure 8 (the radial deformation has been 
enlarged for clearness). 

A 

I J 

B 

C 
D 

Figure 8. Cylinder shapes in different total elongations for the case a = 0.04: (a) 
A = 0.035882224; (b) A = 0.043065200; (c) A = 0.060283469; (d) A = 0.085396683. 

Figure 8(a) represents a low-strain state, Figure 8(b) represents a transition state, 
Figure 8(c) represents a two-phase state, and Figure 8(d) represents a high-strain 
state. As the end displacement increases these four states appears consecutively. 
Here, the asymptotic solutions we obtained can describe the whole deformation 
process. We know that the dimensionless engineering stress 7 represent the stress 
density acting on the cross section of the cylinder in the reference configuration. 
Now, based on the solutions obtained, we can also calculate the true stress density 
a acting on the cross section of the cylinder in the current configuration. Figure 9 
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shows the distributions of the true stress cr corresponding to the four points A — D. 
We see, very interesting, the true stress actually decreases from A to B and to C. 
We also point out that a pure one-dimensional model cannot yield the true stress 
distribution. 
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Figure 9. The distributions of the true stress density a acting on the cross section of the 

cylinder in the current configuration. 

If the diameter-length ratio a is relatively small, we know that there is a snap- 
back on the corresponding 7 — A curve (cf. Figure 6(a) and 6(b)). Thus, for the total 
elongation A located in some special region, there may exist multiple corresponding 
stress values 7, i.e., the solutions are not unique. Then we need to determine which 
solution is the preferred one. 

We take a = 0.02 as an example. From Figure 10, we can see that when 
Ai < A < A2, there are three possible solutions, which are labeled as Ui, U2 and 
U3, respectively. 
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Figure 10. Snap-back on the 7 — A curve for a = 0.02. 
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To choose the preferred solution from Ui, U2 and f/3, we calculate the total potential 
energy values of these solutions. Denote Vt* the total potential energy of the cylinder 
in the displacement-controlled problem and we have 

+ a'iH^VVl + H2V''Vzz))dZ. 

Here we still choose -ffi = 15, = 5. In Figure 11, we plot the differences of the 
total potential energies between solutions U2-, U3 and solution Ui for Ai < A < A2. 
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Figure 11. Differences of the total potential energies between solutions U2, U3 and 
solution Ui for Ai < A < A2 (a = 0.02). 

From Figure 11, we can see that for Ai < A < A*(= 0.038848990), solution Ui 
is the preferred solution and for A* < A < A2, solution f/3 is the preferred solution. 
The engineering stress-strain curve corresponding to the preferred solution is shown 
in Figure 12. 
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Figure 12. The engineering stress-strain curve corresponding to the preferred solution for 

a = 0.02. 

It is worth to note that there exists a jump for the stress 7 at A = A*. This curve 
also captures the main features of the experimental results in the loading process 
(Shaw & Kyriakides 1995, Sun et al. 2000, Tse & Sun 2000, Favier et al. 2001 and 
Li & Sun 2002). 

Figure 13 shows the true stress-strain (c-A) curves (occurring at Z = and 
Z = 0.6) corresponding to the preferred solutions for a = 0.02. 




(a) (6) 

Figure 13. The true stress-strain curves, (a) occurring at Z = for a = 0.02; (b) 

occurring at Z = 0.6 for a = 0.02. 

It is worth to note that for the case occurring at Z = 0, the drop for a become 
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very small and after that there is also a stress plateau (with the true stress value 
larger than the Maxwell stress). For the case occurring at Z = 0.6, the drops for cr 
become larger. But now the stress plateaus have been divided into two parts, which 
corresponding to the fact that phase transition first start at the middle part of the 
cylinder and transform gradually towards the two ends. 

8 Conclusions 

We consider phase transitions induced by axial tension/extension in a slender 
circular cylinder composed of a general compressible hyperelastic material with a 
non-convex strain energy function under clamped end conditions. In order to capture 
the macroscopic phenomena observed in the experiments by others, the problem is 
formulated in a three-dimensional setting and as a result it is governed by a system of 
coupled nonlinear PDE's with complicated nonlinear boundary conditions. A novel 
approach involving coupled series-asymptotic expansions is developed to derive the 
normal form equation (NFE) of the original system. By conducting the phase-plane 
analysis on this NFE, we manage to deduce the global bifurcation properties of its 
solutions. The solutions (including those for post-bifurcation solutions) for both a 
force-controlled problem and a displacement-controlled problem are obtained. These 
asymptotic solutions demonstrate the essential features of phase transitions in a 
cylinder and are consistent with the experimental data qualitatively. Specifically, 
the engineering stress-strain curve is shown to have the same features as observed in 
experiments. These qualitative agreements with experiments give strong supporting 
evidence that a non-convex strain energy function can be used to describe phase 
transitions and the non-convexity of the strain energy function is the main cause of 
the instability phenomena associated with phase transitions in a slender cylinder. An 
important finding is that under the clamped end conditions, there exists a critical 
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stress value 7o(a), which depends on the radius-length ratio and larger than the 
Maxwell stress, such that phase transition can only happen when the engineering 
stress is larger than this value. This is different from the case of the nature boundary 
condition that phase transition can happen as long as the engineering stress is larger 
than the Maxwell stress. 
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Appendix A Incremental elastic moduli 

For initially isotropic material, in the case that there are no prestresses, $ 
should be a function of the principle stretches Ai, A2 and A3, namely $=$(Ai, A2, A3). 
Denote by $j = ^|ai=A2=A3=i, then $1 = $2 = ^3 should vanish since there are no 
prestresses. 

The non-zero first order incremental elastic moduli can be written as 

6 = = "^"11, 

^2 = 0,1122 = *^12, 
^3 = 0-1212 ^ 2^'^^ ~ 
^4 = 0^221 = ^3- 

There are only two independent constants among ^j. 
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The non-zero second order incremental elastic moduli can be written as 





- 

~ "mill 


— y^lll, 






_ 2 

— ^1 111 


= '^'112, 




V3 


2 

— ^^ 1 2233 


= $123, 






2 

— O1II212 


= ^(26^ 


- 2^3 + r]i- 772), 


V5 


2 

— %31212 


4(6 + 


^2-^3), 


V6 


2 

~ '^121323 




^5), 


V7 


2 

— "^111221 


= r]4-^2 




V8 


2 

~ %31221 


= V5- 6, 


V9 


2 

~ "^123123 


= % - C3 





There are only three additional independent constants among rji. 

The non-zero third order incremental elastic moduli can be written as 

$1111, 
$1112, 
•^•1122, 

•^•1123; 

-^(66 + 6^3 - Svi - 3r/2 - 2^1 + 2^2), 

-^(66 + 66 - 6r/2 -Oi- 202 + 8^3), 
-^(6-^2-^/3-^2 + ^4), 
^(66 + 66 + 6Vi - 6r72 + ^1 - 4^2 + S^s), 



61 — aiiiiiiii — 

n 3 

C2 — "11111122 — 

^3 = C^11112222 — 

^4 = C^11112233 — 

^5 = '^11111212 — 

^6 = '^11112323 — 

^7 = '^11221212 — 

^8 = ^^11221313 — 

^9 = '^12121212 — 
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^10 - a?2121313 - q^9' 



6 

On = a?n2i323 = "^(^^i " + 3r/3 -9^ + 6, + 39, - 36,), 
Oi2 = a?nii22i = ^(66 + 66 - 3r/i - 3r/2 + 2^1 - 2^2), 

6*13 = «11112332 = -2(^72 - ^3 + ^4), 

^14 = a?i22i22i = ^(66 + 6^3 - 6r/2 + + 2^^2 - 3^3), 

Ol5 = «11221331 = ^(6 - ^72 - ^73 + ^2 - O4), 

Om = a?2i2i22i = -^(66 + 6^3 - + 4^2 - 3^3), 

o 

On = a?22ii22i = ^(66 + 66 - 2r/i + 2r/2 + ^1 - 4^2 + 3^3), 

6*18 = ^12121331 = 3^16, 

^19 = ai22ii33i = ^(6^1 ~ 3?7i - 3?72 + 6?73 + 6*1 - 46*2 + 3613), 
^20 = ani23i23 = ^(66 + 36 - 3r/i + 3r/3 + ^1 - ^2 - 3^3 + 3^4), 
^'21 = ani23i32 = ^20 - ^(6 + 6 - ^1 + ^2), 
^22 = a?2i22323 = ^(66 - 66 + 12r/2 - 12r/3 + ^1 - 4^2 + 3^3)- 
There are only four additional independent constants among 9i. 

The expressions of the constants q (i = 1, ■ ■ ■ , 12) in equation (5.3) are given 
below: 

6^ 6^ ^|r/2 , ^ 6^ , 6^ ^1^3 _ ^ _ 56^4 
4 46 46 466 4 46 46 466 4 46 
_ 3^|r/4 _ 5^fr/4 _ 36^4 _ HIva ,V5_ 6^ _ ^j^s 6^ 

86' 866' 46 466 4 46 86^ 866' 46 



466 4 46 46 466 4 46 46 466 
_ 6^2 _ 6j/2 6^ ^1^2 6^ _ 6% ^i^3 , !/i _ 56^4 66^4 
46 46 26 466 46 46 466 2 46 46^ 

_ ^ _ 5^fr/4 6^ _ 5^|r?4 _ 6^ 66^5 _ ^f^s 6^/5 _ ^1^5 _ ^ 

46' 866' 46 466 46 46^ 866' 46 466 2 

^ 56^7 _ 6^ ^ 5^ir/7 ^ 6^ _ 6^ ^ ^|r/8 



46 46 466 46 46 466 
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4^1 866 4^26 4 46 26' 46 866 46^6 

4 26 46 266 2 26 4^1 466' 266 2 
26 4^1 466' 266 2 26 266 2 26 266' 

_ _ ^IVl ^2V2 _ §J]2 _ 6^/2 _ 6^ ^2^ _ ^lV2 ^2Jh _ ^lV3 

26' 46'6 ^1 ei 46 26 66 26'6 6 46 

3^2^3 _ ^ _ 6^ , 66^4 ^1^/4 _ ^iV4 6^ _ ^ _ 6^ 66^5 

46^3 2 6 4e| 46' 266' 46 6^3 26 46' 
_ ^1^5 6^/5 _ ^1^5 , !/7 6^ _ 6j/7 6^ _ ^im ^'^s 

46e| 46 266 2 + 6 46 6^3 26 46 266' 

_ ^ _ 6^ 36^2 _ ^r]2 _ 6^ 5^|r/2 _ ^2V2 _ 6^ 66^4 

8 86 46 26' 46 866 46'6 26 46' 
_ ^2^4 6^ _ ^im 6^ _ 6j/7 ^1^7 
46el 46 266 26 46 266' 

?72 6^ 6^ ^1^2 6^ _ !^ _ _ 66^4 _ ^|?74 _ ^iv4 

2 26 46 266 46 2 26 26' 46' 466' 
_ 6r/4 _ 6^ _ ^iv4 ,vi, 6^/7 6^ 6^ ^'^7 

26 26 266 2 26 26 26 266' 
6^1 6^ _ 6^ Cl^2 _ 6^ ^ _ ^ _ ^iv4 6^ _ 6^ 
46 26 46 266 26 86' 26' 466' 46 26 
_ ^2^4 6^ _ 6^ 6^ 

266 86' 26 46 26 266' 

Hi + _ + ^2^1 , 6^/2 _ 56'^4 66^4 _ Cf^ _ 6^ _ 6^ 
8 86 86 866 46 326' 326' 86' 86 86 
_ 66^5 6^ 6^ 

326' 326' 86 86' 

6^1 ^i^i _ 366^4 _ ^r]4 _ 6^ 66^5 6^ 

86 866 166' 326' 86' 86 326' 86' 

_m_ 6^ _ 6^ _ ^2Vi _ I ^'^2 6^ 66^4 

16 86 166 866 16 86 166 866 86 86' 
_ ^|r/4 _ 6^ _ 6^ 66^5 _ _ 6^/5 _ 6^ 6^ 6^ 

86' 86 86 86' 86' 86 86 86 86 

6^ 6^ 

86 86 ' 

Hi _ _ ^^1^2 6^ + !^ + + ii!^ + 66^4 _ _ 

8 86 86 46 8 46 86' 86' 46' 86 

_ 6^ _ f^fe 6^ 6^ 6^ 6^ 
46 86' 86' 86 86 46 86 ' 
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86 166 866 16 86 166 866 8^1 8e| 86 86' 
Appendix B Non-dimensional field equations 

The full forms of the non-dimensional field equations (3.3) and (3.4) are given 
below: 

(26 + 26)t^. + ^iw,, + 46w. + s((26 + 26)?;.. + ^^sWss) 
+ e[(2?72 + 2?73 + 2r]7 + 2r]s)vv^ + (2?72 + 2r]7)v^Wz + 2r]2VWzz + Vi^zW^z 
+ (4?74 + 4?75)tw, + Arj^WzWs + s{r]iVzVzz + (2?72 + Sr/s + 107/7 + 2r]s)vsVz 
+ 2r]2WzzVs + 2ri'jv^^w, + (20?74 + 'ir]5)vsWs + (2?72 + 2?73 + 2?77 + 2r]s)vVsz 
+ {2r]2 + 2r]j)wsVsz + '^r]7V^Wsz + ^V^WsVOsz + (4?74 + 4?75)f ty.. + 4?74ty^tys. 
+ s^((4?72 + A'q^)vsVsz + ^rijVzV^s + 8?74w;st'ss + 8?74WsW;ss)] + e^[(6'2 + 3614 + 6^12 

+ + 2^15)1;'?;. + (2^3 + 2^4 + 2^14 + 2d^^)vVzW, + {02 + d^2)VzWl + (^3 + ^4)^'^;.. 

+ 2e2VWzWzz + ^^iw;>.. + (2^5 + 2^6 + 408)?;^u;s + (4^7 + ^e^)vWzWs + 205U;>, 
2 11 1 

+ ■5((-6'i6 + -6*7 + -6'8)?;f + {9^ + 6^)VV^V^^ + O^V^WzV^z + ^O^V^Wzz 

+ (10012 + 2013 + 12015 + 202 + 604)t;t;.t;. + (lO0i4 + 20i5 + 203 + 29i)vzWzVs 

+ (203 + 2e4)vWzzVs + 2e2WzWzzVs + (2014 + 2015 + 40i7)t;>s + (2014 + 20i5)w,,W, 

+ 2012^^^..^, + 20i2t;.w^.M;, + (2005 + 406 + 2408)t;t;,'«;, + (2007 + ^Os)wzVsWs 

16 

+ (8016 + 207 + 29s)VzW^, + 295WzzWl + —9gwl + (012 + 013 + 2015 + O2 + Wa)v\sz 

+ (2014 + 2015 + 203 + 2ei)vWzVsz + (012 + 02)wlvsz + (4014 + 40i5)i;t;,u;,, + 2e^wlwss 

+ 4:9i2VzWzWsz + (807 + 89s)vWsWsz + Se^WzWsWsz + (205 + 206 + 4:9s)v'^Wss 

+ (407 + 4es)vWzWss) + s\2erVzVzzVs + (16012 + 4015 + 202 + 2e4)v,vl + IQOrVsWsWsz 

+ 2e^WzzVl + ^euVzzVsWs + (016 + Oj)vlVsz + (4012 + 4015 + 402 + ^9i)vV,Vsz 
+ (4014 + ^0^)WzVsVsz + (4014 + ^0^7)VzWsVsz + (4016 + 407)?^^'.. + ^O^iVzVsWsz 



42 



ss 



+ (4:012 + 4:6i5)vVzVss + 4:6uVzW^Vss + {865 + 86s)vWsVss + SdjWzWsVss + 26i7vlw 
+ {865 + 8ds)vVsWss + 867WzVsWss + 8diQV^WsWss + 86gw1wss + (326^5 + 8ds)vlws) 
+ S%{4ei2 + 4e2)vlvsz + S^isf^Ws^^.s + lQe5VsWsV,, + 8e5vlwss)] = 0, 

(SI) 

5 1 

^3Vzz + 8^iVs + (2^2 + 2^3)^^^ + 4:S^iVss + ^[(2^4 + -^VbWz + iv4 + Vb)vv^^ 

+ V^WzV^z + V^v^w^z + (8771 + 8r]2)vvs + 8r]2WzVs + 8ri7V^Ws + 2r]7WzzWs 

+ (6?74 - 2r]5)wl + (2?72 + 2r]3 + 2777 + 2r]^)vWsz + (2?72 + 2r]7)wzWsz + s{2r]4VzzVs 

+ (14?7i + 2?72)f^ + irj^v^Vsz + 8r]iWsVsz + (4?72 + A7]^)vsWsz + (4?7i + 4?72)fWss 

5 1 

+ Ari2W^Vss + 4r]'jV^Wss + 8r]4WsWss) + 8s^?7iWsWss] + ^^[(g^s + + 36*8)^^^ 
+ (^7 + ^8)f f^^f.. + 0^v,w,w,, + (4^1 + 8^2 + 4e^Wv, + (8^2 + 8ei)vw,Vs 

VVzWs + 86iiVz'Wz'Ws + (26'i4 + 26i^)vWzzWs + 26i2WzWzzWs 
+ (6^5 - 2^6 + 4^8)t^M^' + (6^7 - 2es)w,wl + (^12 + ^13 + 2^i5 + ^2 + 3^4)^^'^,, 
+ (2^14 + 2^15 + 2^3 + 2e^)vW,Wsz + (^12 + 02)wlwsz + s(^^9^>.. + (8^5 + 2^8)v> 
+ (2^5 + 2es)vv,,Vs + 2^7w;^t;,,t;, + 2e7V,w,,Vs + (14^i + 16^2 + 2^3)^,^ 
+ (146'2 + 2ei)w^v'l + 2^i6f^W2^Ws + (286'i2 + Aei^)v^VsWs + ^OuW^^VsVOs 
+ 26'i7f^^w^ + 2Ae^Vsw'l + (46'5 + 46'8)w^^;s^ + 46'7t;^w^t;s^ + (86'i2 + 8di^)vWsVsz 

+ 8^i4tf^M;,i;,^ + (^16 + e'j)vlWsz + (4^12 + 4^i5 + 4^2 + 4^4)^^^^sU^s^ 

+ (4^14 + Ae^)w,vsw,, + (4^14 + Aen)v,WsWs, + (4^16 + 4ej)wlw,, 

+ (2^1 + 4^2 + 2^3)^^^^ss + (4^2 + 4^4)^^^^^^,, + 2e^wlvss + (4^12 + Aei^)vv,Wss 

+ A9iaVzW:,Wss + (86*5 + 89s)vWsWss + 89iWzWsWss) + s^{20^Vzzvl + (y6'i + ^612)^^ 

+ 86'5f^t;^t;^^ + 166'i2t'sU'st;52 + (46*12 + 402)vlwsz + 2e^vlvss + (86*1 + 86'2)wWsWss 

+ 85=^^1^72^;,,] = 0. 

m 
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Appendix C 



The full forms of the terms Hi-Hj in equations (4.4)-(4.8) are given below: 

H, ={26, + 2^6 + ^es)Vo^Wi + {Ou + ^13 + 2^15 + ^2 + 3^4)K)'K). 

+ (4^7 + 4es)VoWiWo, + (2^14 + 2^15 + 2^3 + 2e,)VoVo.Wo, + 2e,WiWl 

+ {612 + e2)VoM, + (^3 + ei)Vo^Wo.. + 2e2VoWo.Wo^, + \bxWIw^,,, 

i72 =(24^5 + 8^6 + 3208)^0^1^^1 + -ttQ'aWI + (8^5 + 8^6 + m^)V^W2 

o 

+ (12012 + 4013 + 16015 + 402 + 1204)\/o\/i\/o, + (8016 + 207 + 208)W^,Vo, 
+ (2014 + 2015 + 40i7)W^in', + (^016 + \Qi + \Q^)yl + 805^^21^0'. 

+ (2012 + 2013 + 4015 + 202 + ^B^V^Vxz + (2407 + 808)V^iiyil^o. 

+ (1607 + 1608)nW^2W^o. + (12014 + 4015 + 403 + 404)V^iV^o.Vro, + (4014 + 4015 

+ 403 + 404)K)V^i.Vro. + (12012 + 202)1^1.1^0'. + (12^7 + 1208)l/oW^iW^i. 

+ (6014 + 6015 + 203 + 204)K)K).W^1. + 1205^^11^0.1^1. + (6012 + 202)^.1^0.1^1. 
+ (2014 + 20i5)V^oVri^.. + (07 + ^8)K)14,.K).. + 201211^1 W^o.V^O.. + ^5K).11^0.V^0.. 

+ (403 + 404)1^0^11^0.. + 205l^i'Vro.. + 201211^1 l^o.VTo.. + ^^5^' H^o.. 

+ 40214 11^0. 11^0.. + 202l/oll^l.ll^O.. + 01 11^0. 11^1.1^0.. + (^3 + 04)l^o'W^l- 
+ 202l^oll^o.ll"i.. + ^^1 11^0.^1-, 
=(401 + 802 + 403)^0'^! + (6^5 - 206 + 408)K)ll^i' + (8012 + 80i5)ll^il/o. 

+ (^05 + ^06 + 308)l/oK)' + (802 + 804)l/ol411^O. + (607 " 208)ll^i'iyo. 
+ 80141^11^0.11^0. + ("-h + \h)ylW^z + 4031^1 IVo', + (012 + 013 

+ 2015 + ^2 + 304)Vo'W^i. + (2^14 + 2015 + 203 + 204)l/oVro.l^i. 

+ (012 + 02)11^0.^1. + (i^5 + ^^6 + 08)K)'^O.. + (^7 + 08 ) 1^0 W^O. H.. 
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Ha ={\oi + + ^^3)v;,' + {\e2 + ^^4)i^o'w^o. + (^^3 + \eA)v,wl + ^e^w^o'., 

=(3r/i + 5r/2)^^i + '^rfS^I + 27/7l^iVo. + ^^4^)' + (3^2 + V^)ViW^z 
+ (^2 + r/3)VoVri2 + r/2Vro^iyi^, 

+ ^6 + 2^8)^o'W^l + {\0l2 + ^^13 + ^15)K)'^0. + (2^7 + 2^8)^0 W^lW^O. 

H-j =(6r/4 + 2r/5)^iW^i + (4r/4 + 4/75)^^^2 + (3r/7 + r/s)^!^. + {m + r/s)^^^!. 
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